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1 .  INTRODUCTION 


Ford  and  Fulkerson  (1956)  observed  that  an  easy  way  to  find  the 

value  of  the  maximal  possible  flow  through  an  undirected  source-sink 
* 

planar  network  is  to  construct  the  dual  graph,  assign  the  capacities 
of  the  intersected  primal  arcs  as  lengths  of  the  undirected  arcs  in 
the  dual,  and  then  find  the  shortest  route  through  the  dual.  The 
length  of  this  shortest  route  will  be  equal  to  the  value  of  the  mini¬ 
mum  cut  set  of  the  primal  network  and  hence  is  equal  to  the  value  of 
the  maximal  flow  because  of  the  min-cut  max-flow  theorem. 

Extension  of  this  idea  to  directed  networks  began  with  de 
Ghellinck  (1961) .  He  developed  a  convention  for  constructing  a  dual 
of  a  directed  source-sink  planar  network  but  was  not  concerned  with 
arc  lengths.  In  a  paper  written  in  1969,  Doulliez  and  Rao  (1971) 
provide  the  dual  construction  and  the  arc  length  conventions  for 

I 

directed  capacitated  flow  networks  when  zero  lower  bounds  exist  on 
arc  flows.  This  convention  was  apparently  the  result  of  conversations 
with  de  Ghellinck. 

Because  of  a  need  for  dual  graphs  of  directed  network  in  studies 
of  network  interdiction  models,  the  author  (1970)  independently  devel¬ 
oped  tf  j  same  convention  as  Doulliez  and  Rao  but  went  one  step  further 
and  provided  a  convention  for  handling  problems  where  the  lower  bound 
on  arc  flow  could  take  any  negative  value  or  any  positive  value  not 
exceeding  the  upper  bound  capacity.  This  paper  also  observed  that  for 

Fold  and  Fulkerson  (1956)  defined  source-sink  planar  networks  as  planar 
networks  which  remain  planar  after  an  arc  connecting  the  source  and  sink 
is  added  to  the  network. 

/- 


2 


flow  networks  having  positive  lower  bounds  on  arc  flows,  an  infeasible 
primal  problem  is  indicated  by  a  cycle  of  negative  length  in  the  dual 
network.  In  a  later  paper  (1971)  the  author  proved  that  a  cycle  of 
negative  length  was  a  necessary  and  sufficient  condition  for  infeasi¬ 
bility.  In  addition,  a  proof  was  provided  that  the  value  of  the  minimal 
feasible  flow  could  be  obtained  by  changing  the  sign  on  the  arc  length 
of  the  dual  and  finding  the  longest  route  from  the  destination  back 
to  the  origin. 

During  1970,  M.  F.  Sakarovitch  was  also  working  on  these  same 
problems.  His  research  (1970)  was  motivated  by  some  conversations 
with  de  Ghellinck.  He  also  proved  that:  a  cycle  of  negative  length  was 
a  necessary  and  sufficient  condition  for  infeasibility.  In  addition, 
he  proved  that  an  optimal  flow  in  a  primal  arc  of  a  maximal  flow  net¬ 
work  is  equal  to  the  difference  between  the  dual  shortest  distances 
from  the  origin  to  the  two  nodes  at  opposite  ends  of  the  dual  arcs 
intersecting  the  primal  arc.  This  important  property  will  play  a 
major  role  in  the  algorithm  to  be  presented  below  for  non-planar 
directed  networks. 

All  of  the  work  just  summarized  is  restricted  to  those  networks 
which  are  source-sink  planar.  The  main  purpose  of  this  paper  is  to 
extend  these  ideas  to  general  flow  networks. 

In  the  sections  to  follow,  we  will  first  review  the  "state  of 
the  art"  for  source-sink  planar  directed  network  flow  problems.  We 
will  then  turn  to  the  problem  of  non-planar  networks.  A  pseudo  dual 
will  b'i  defined  and  tha  special  form  of  the  dual  shortest  route  problem 
will  be  stated.  The  details  of  the  algorithm  for  solving  this  problem 
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will  then  be  presented  followed  by  an  example,  and  the  proofs  which 
validate  the  algorithm.  The  final  section  presents  a  brief  discussion 
of  the  usefulness  of  the  algorithm. 


- Uk — 
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2.  MAXIMAL  FLOW  IN  A  PLANAR  NETWORK 

The  original  flow  network  will  be  called  the  primal  network. 

A  mesh  of  a  planar  primal  network  is  any  region  surrounded  by  nodes 
and  arcs  but  containing  neither  in  the  plane  on  which  the  network  is 
constructed.  The  region  of  the  plane  completely  surrounding  the 
primal  network  will  be  called  the  external  mesh. 

The  construction  of  the  dual  of  a  source-sink  planar  directed 
network  consists  of  the  following  steps  (McMasters  (1970)): 

1.  Connect  an  artificial  arc  between  the  sink  and  source 
of  the  primal  and  position  it  below  the  network  such  that  it 
crosses  no  other  arcs.  The  resulting  network  will  be  referred 
to  as  the  modified  primal  network. 

2.  Place  a  node  in  each  mesh  of  the  modified  primal  in¬ 
cluding  the  external  mesh.  Let  the  origin  of  the  dual  be  the 
node  in  the  mesh  involving  the  artificial  arc  and  the  destina¬ 
tion  be  the  node  in  the  external  mesh. 

3.  For  each  arc  in  the  primal  (except  the  artificial  arc) 
construct  two  oppositely  directed  arcs  that  intersects  it  and 
join  with  nodes  in  the  meshes  adjacent  to  it. 

4.  Assign  the  value  of  the  upper  bound  capacity  of  the 
primal  arc  as  the  length  of  the  intersecting  dual  arc  having 
the  same  direction  that  the  primal  arc  would  have  if  it  were 
rotated  90  degrees  counterclockwise.  Assign  to  the  oppositely 
directed  dual  arc  a  length  equal  to  the  negative  of  the  lower 
bound  capacity  of  the  primal  arc. 
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The  convention  described  in  step  4  is  illustrated  in  Figure  1 
where  the  arc  (i,j)  is  the  primal  arc  having  its  flow  constrained 

between  two  values  and  such  that  .  The 

adjacent  dual  nodes  are  u  and  v.  Arc  (u,v)  is  assigned  a  length 


equal  to  while  arc  (v,u)  is  assigned  a  length  equal  to  -L^ . 

We  assume  M^.  s  0  but  place  no  restriction  on  L^_.  except  that  it 
cannot  exceed  . 


Figure  2a  is  an  example  of  a  flow  network  having  node  1  as  its  source 
and  node  4  as  its  sink.  The  numbers  on  each  are  represented  , 

.  Figure  2b  shows  the  initial  phase  of  the  construction  of  the 
dual.  The  dual  nodes  are  designated  as  A,  B,  C  and  D.  Each 
dashed  arc  corresponds  to  the  location  of  each  pair  of  dual  arcs.  Fig¬ 
ure  2c  is  the  completed  dual  with  its  associated  arc  lengths. 


Ford  and  Fulkerson  (1962)  defined  the  generalized  value  of  a 
primal  cut  set  when  lower  bounds  are  positive  by  equation  (1)  where 
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the  set  S  contains  the  primal  source  node  and  the  set  S  contains 
the  primal  sink  node. 


V(S,S) 


l 

ieS 

jeS 


M. 


l  L.. 

icS  1J 
JeS 


Oj 


They  then  extended  the  min-cut  max-flow  theorem  to  include 
this  definition  of  the  value  of  a  cut  set;  namely,  if  a  feasible  flow 
can  be  found  for  a  network  with  nonnegative  lower  bounds  on  arc  flows 


then  the  maximal  possible  flow  through  the  network  is  equal  to  the 
minimal  generalized  value  of  all  possible  cut  sets  of  the  (S,S)  type. 
Later  Fulkerson  (1962)  dropped  the  restriction  that  the  lower  bound 
be  nonnegative. 

The  relationship  between  the  convention  of  Figure  1  and  the 
generalized  value  of  a  primal  cut  set  is  described  by  the  following 
theorem: 

The.oA.ejv  1 :  Any  dual  chain  without  cyclei  directed  ^ rom 
the  dual  source  to  the  dual  -ilnk  ha6  a  length  equal  to 
the  generalized  value  o&  the  primal  cut  i>ct  Intersected 
by  the  arc 6  o<$  the  chal.n. 

Proof:  Each  path  without  circuits  through  the  dual  of  a  source- 
sink  planar  undirected  graph  intersects  a  cut  set  of  the  primal  network 
such  that  the  primal  source  and  sink  are  separated  by  the  cut  and 
there  is  a  separate  path  corresponding  to  each  cut  (Ford  and  Fulkerson 
(1956),  Whitney  (1933)).  This  is  still  true  if  the  primal  and  dual 
arcs  are  replaced  by  directed  arcs.  Further,  because  tho  dual  graph 
under  the  Figure  1  convention  has  two  oppositely  directed  arcs  between 
every  pair  of  nodes  in  adjacent  primal  meshes  it  is  also  true  that  a 
separate  chain  exists  directed  from  the  dual  source  to  the  dual  sink 
corresponding  to  each  primal  cut. 

Suppose  we  select  an  arbitrary  cut  set  of  the  primal  network 


such  that  the  source  is  in  a  set  S  and  the  sink  is  in  set  S.  Any 
such  cut  set  can  be  illustrated  as  shown  in  Figure  3  for  planar  graphs. 


Notice  that  this  set  contains  both  arcs  directed  from  S  to  S  and 
from  S  to  S.  Let  the  chain  directed  from  the  dual  source  to  the 
dual  sink  which  intersects  the  cut  set  arcs  be  a-b-c-d-e-f. 


Figure  3. 

Using  the  convention  of  Figure  1  we  assign  a  length  of  -Lfc 
to  dual  arc  (a,b)  which  crosses  primal  arc  t,  a  length  of  Mg  to 
arc  (b,c)  which  crosses  arc  s,  a  length  of  to  arc  (c,d), 

a  length  of  -L^  to  arc  (d,e),  and  a  length  of  to  arc  (e,f). 

The  chain  directed  from  a  to  f  would  therefore  have  a  length  of 


M  +  M  +  M  -L  -  L 
s  r  p  t  q 

which  we  see  immediately  corresponds  to  the  form  of  the  right-hand 
side  of  (1)  for  this  arbitrary  cut  set  (the  arcs  p,  r  and  s  have 
ieS  and  jeS  while  arcs  q  and  t  have  ieS  and  jeS).  Because 
our  cut  was  arbitrary  we  realize  that  if  we  analyze  any  cut  set  of  the 
(S,S)  type  in  this  manner  we  will  get  the  length  for  its  dual  chain 


being  described  by  this  general  form.  Therefore,  each  chain  without 
cycles  from  the  dual  source  to  the  dual  sink  will  have  a  length  equal 
to  the  generalized  value  of  the  primal  cut  set  it  intersects. 

As  a  consequence  of  Theorem  1  and  the  generalization  of  the 
min-cut  max-flow  theorem  we  get  the  following  corollary: 

CoKoWx/iy  1:  The  value  o maximal  filow  in  a  fieotible 
diAected  AouAce-Aink  planaA  neiwoAk  u>  equal  to  the 
length  otf  the  ihoateAt  aouie  through  the  dual  neMvoak. 

Recalling  Figure  2a  we  know  that  the  value  of  the  primal  cut 
set  consisting  of  arcs  (1,2),  (2,3),  and  (3,4)  is 

M12  +  M34  -  l23  =  5  +  3  *  1  =  ?•  (2) 

The  route  through  the  dual  network  which  corresponds  to  this  cut  set 
is  the  chain  [  (A,C) ,  (C,B) ,  (B,D)  ]  which  has  a  length  given  exactly 
by  equation  (2) .  This  chain  happens  to  be  a  minimum  route  of  Figure 
2c  so  that  the  maximum  possible  flow  through  the  primal  network  is, 
in  fact,  7  units.  As  we  will  show  below,  the  fact  that  we  found  a 
shortest  route  of  finite  length  means  that  the  primal  problem  is  feasible. 

The  extension  of  the  min-cut  max-flow  theorem  involving  the  use 
of  the  generalized  value  of  a  cut  set  as  given  by  (1)  is  contingent 
on  the  existence  of  a  feasible  flow.  Feasibility  requires  that  a  set 
of  arc  flows  can  be  found  which  satisfy  the  flow  conservation  equations 
at  each  node  as  well  as  the  upper  and  lower  bounds  on  individual  arc 
flows.  When  the  lower  bound  is  nonpositive  then  we  know  that  a  flow 
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problem  will  always  be  feasible  since  X 


ij 


0  is  a  feasible  solution 


under  the  ^  0  assumption.  If,  however,  the  lower  bounds  on  some 

arcs  are  positive  then  we  cannot  claim  such  problems  will  be  always 
feasible.  An  example  of  an  infeasible  problem  can  be  obtained  by 
changing  the  lower  bound  on  arc  (1,3)  in  Figure  2a  to  a  value  of 
3. 

The  answer  to  the  feasibility  question  is  provided  by  Theorem  2. 


Theorem  2:  A  necehhojiy  and  hu^ident  condition  {^on.  an 
infieaiible  £ low  problem  a  ho  dated  with  a  dixected 
houAce-hink  planaA  network  ih  that  a  cycle  o negative 
length  exiht 5  in  the  dual  network. 


The  proof  of  Theorem  2  depends  upon  the  circulation  theorem  of 
Hoffman  (1960)  which  is  addressed  to  networks  having  a  sourceless- 
sinkless  structure.  Such  a  network  is  easily  created  from  any  flow 
network  by  connecting  a  directed  arc  from  the  sink  to  the  source  and 
specifying  the  appropriate  upper  and  lower  bounds  on  its  arc  flow. 

The  circulation  theorem  states  that  a  network  flow  problem  will 
be  feasible  if  and  only  if 


I 

icY 

jeY 


it 


L  L 

ieY 

jeY 


ij 


(3) 


for  all  cut  sets  (Y,Y)  of  the  network  after  it  has  been  converted 


to  the  sourceless-sinkless  form. 
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Suppose  we  define  the  value  of  the  cut  set  (Y,Y)  as 

V(Y,Y)  -  l  M  -  £_  L  .  (A) 

ieY  J  ieY  J 
j  eY  jeY 

If  we  do  this  then  the  feasibility  condition  given  by  (3)  is  equiva¬ 
lent  to  V(Y,Y)  ^  0  and  infeasibility  corresponds  to  some  V(Y,Y)  <  0. 

An  interesting  correspondence  between  dual  cycles  and  V(Y,Y) 
is  described  by  Theorem  3.  We  will  prove  this  theorem  before  proving 
Theorem  2  since  it  will  be  useful  in  the  proof  of  the  latter. 

Theorem  3:  JheAe  l&  a  ^epaicrfe  dual  cycle,  intent,  ectlng 
(Uict>  oft  each  pAlmat  cut  &et  (Y,Y).  Thli  cycle  haA  a 
length  equal,  to  V(Y,Y). 

Proof:  Any  circuit  in  the  dual  corresponds  to  a  cut  set  (Y,Y) 
of  the  primal  (Whitney  (1933)).  Because  of  our  convention  of  having 
two  oppositely  directed  dual  arcs  intersect  each  primal  arc  we  know 
that  two  cycles  exist  (oppositely  directed)  in  the  dual  corresponding 
to  each  cut  set  of  the  (Y,Y)  type  in  the  primal.  Suppose  the  set 
Y  is  circumscribed  by  these  dual  cycles.  The  convention  of  Figure  1 
will  then  result  in  the  cycle  directed  clockwise  having  a  length 
corresponding  to  V(Y,Y).  If  Y  is  the  set  circumscribed  by  the 
cycles  then  the  counterclockwise  cycle  has  a  length  V(Y,Y). 

Figure  4  is  an  example  of  an  arbitrary  cut  set  of  the  (Y,Y) 
type  with  its  associated  dual  cycles.  Primal  arcs  e,  f,  h  and  j 
are  directed  from  Y  to  Y;  arcs  g,  i,  k  and  l  are  directed 
from  Y  to  Y.  The  convention  of  Figure  1  will  result  in  the  cycle 
p-q-r-s-t-u-v-w-p  having  a  length  of 


L2 


M  +  M.  +  M.  +  M  -L  -  L,  -  L,  -  L. 
g  i  Tc  e  f  h  j 


which  is  the  right  hand  side  of  (A) . 


We  will  now  prove  Theorem  2.  Suppose  that  the  primal  problem 
is  infeasible.  From  the  circulation  theorem  we  know  that  there  will 
be  some  cut  set  (Y,Y)  of  the  sourceless-sinkless  primal  which  has 


l  Mij  <  L  L 


ieY 


ieY 


ij' 


(5) 


jeY 


jeY 


We  can  rewrite  (5)  in  the  following  form: 


l  M..  -  l  L,,  <  0. 
ieY  ij  ieY  ij 
jeY  jeY 


(6) 


from  (A)  we  know  that  the  left-hand  side  of  (6)  is  the  value  of  the 
cut  set  (Y,Y).  It  is  also  the  length  of  a  cycle  in  the  dual  because 
of  Theorem  3.  Therefore,  an  infeasible  primal  problem  results  in  a 
cycle  of  negative  length  in  the  dual. 
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To  show  sufficiency,  we  consider  a  cycle  of  negative  length  in 
the  dual.  Let  its  length  be  represented  by  D.  Because  of  Theorem  3, 


»  “  l  M  -  l  L  *  V(Y,Y), 
ieY  J  ieY  J 


and  D  being  negative  means  that  the  value  of  some  cut  set  (Y,Y)  in 
the  primal  is  negative.  For  a  cut  set  (Y,Y)  to  have  a  negative  value 
(5)  must  be  satisfied  and  the  primal  will  be  infeasible. 
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3.  MINIMAL  FLOW  IN  A  PLANAR  NETWORK 


It  is  sometimes  useful  to  also  know  the  minimal  feasible  flow 
for  a  network  having  positive  lower  bounds  on  arc  flows.  The  following 
theorem  provides  an  easy  way  of  determining  such  a  flow. 

Th.exiA.ejm  4:  l&  the,  oac  length*  the  dual  o&  a  fieaAlble 
pJupial  filou)  netwoAk  have  tiieoi  *ign*  AeveA*  ed  then  the 
longest  chain  dinected  fi-iom  the  dual  deAtinatlon  back 
to  the  dual  AouAce  mil  have  a  length  equal  to  the  value 
o&  the  minimum  fieaAible  filou)  thAough  the  netuioAk. 


The  proof  of  this  theorem  is  based  on  the  min-flow  equivalent 
of  the  max-flow  min-cut  theorem  (Ford  and  Fulkerson  (1962));  that  is, 


min  Q  =  max  L  -  ^  M  ,  (7) 

ieS  ieS  J 

jeS  jeS 

where  Q  represents  the  total  flow  leaving  the  network. 

The  combination  of  our  convention  of  allocating  dual  arc  lengths 
and  the  sign  reversal  just  described  will  cause  each  chain  from  the 
destination  back  to  the  origin  in  the  dual  network  to  have  a  length 
equal  to 


l  -  L  Mij  ■ 

ieS  J  ieS  3 
jeS  jeS 

The  proof  of  this  property  parallels  that  of  Theorem  1.  Therefore, 
we  get  the  right-hand  side  of  (7)  when  we  find  the  longest  chain  back 


through  the  dual. 
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A.  OPTIMAL  ARC  FLOWS  IN  PLANAR  NETWORKS 

The  final  aspect  of  the  source-sink  planar  network  analysis  is 
the  determination  of  the  optimal  arcs  flows  of  the  primal  network. 

These  flows  are  available  from  the  solution  to  the  dual  network  shortest 
route  or  longest  route  problem  depending  on  whether  we  want  maximal  or 
minimal  flow  through  the  network. 

Consider  the  maximal  flow  problem.  Let  A  be  the  set  of  all 
dual  arcs  associated  with  a  given  primal  flow  network.  Then  the  dual 
network  shortest  route  problem  can  be  stated  as  follows: 

minimize  £  ^  U±j  , 

f  -1  (i  =  1), 

N 

subject  to  l  [U  -  U  ]  -<  0  (i  =  2,3 . N-l) ,  (8) 

j=l  J  3 

\  +1  (i  =  N), 

and  U±j  2:  0  ((i,  j)eA)  . 

If  we  think  of  U^_.  and  as  representing  respectively  a 

flow  and  a  unit  flow  cost  through  a  dual  arc  (i,j)  then  we  realize 
that  (8)  corresponds  to  a  minimal  cost  network  flow  problem  for  a 
total  flow  out  of  one  unit.  The  optimal  solution  will  yield  a  chain 
flow  of  one  unit  over  the  cheapest  route  from  the  origin  to  the  desti¬ 
nation  of  the  dual  network.  This  route  obviously  corresponds  to  the 
shortest  route  through  the  dual  network.  The  dual  linear  programming 
problem  associated  with  (8)  can  be  written  as 
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maximize  V.,  -  V, , 

N  1 

(9) 

subject  to  Vj  -  £  £  ( (i,j ) cA)  . 

If  we  apply  any  of  the  current  algorithms  (see  Dreyfus  (1969)) 
for  determining  the  shortest  routes  from  node  1  to  every  other  node 
we  will  obtain  a  number  associated  with  the  length  of  the  shortest 

route  from  node  1  tc  node  i.  If  we  let  ■  f^  then  we  will  have 
an  optimal  solution  to  (9)  • 

The  relationship  between  the.  values  and  the  optimal  arc 

flows  in  the  original  primal  network  is  described  by  Theorem  5  (Sakaro- 
vitch  (1970)). 


Theoxem  5:  Confide. A  a.  pnmal  oac  (a,B)  duAected  in.om 

node,  a  to  node  g  ivkich  hcu  a  ilou)  X  .  Let  the 

op 

intersecting  dual  oaca  be  (i,j)  and  (j,i)  uuth  the 
length  o&  (i,j)  being  M  and  the  length  o 6  (j,i) 

ap 

being  ~L  ..  I  ^  X  a  Aati&iie& 

&p  Ct  p 


(10) 


ion.  eveAy  oac  (a,g)  then  an  optimal  fieaAible  solution 
hoj>  been  fiound  to  the  pnlmal  maximal  iloai  pAoblem. 


Proof:  Because  is  a  feasible  solution  to  (9)  we 


know  that 


L  .  S  V  -  V  S  M 
ag  j  i  ag 


under  our  convention  for  assigning  dual  arc  lengths.  Thus 
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L  .  SX  .  iM 
ag  ag  ag 


and  therefore  is  a  feasible  flow  through  arc  (a,p). 

To  check  for  flow  conservation  at  the  primal  nodes  we  first 

introduce  a  new  primal  arc  (n,l)  connecting  the  sink  to  the  source 

and  set  its  upper  and  lower  bounds  as  well  as  the  arc  flow  from  the 

sink  to  the  source  at  the  value  of  V.  -  V,  .  We  then  have  a  source- 

N  1 

less-sinkless  structure  and  the  flow  conservation  equations  can  be 
written  as 

(ID 


y  fx  „  -  x„  i  =o 

“  L  <*3  6a 


3 


for  ot  1,2, .. .  ,n. 

If  we  next  add  dual  arcs  to  intersect  the  primal  (n,l) 
arc  we  will  create  a  dual  network  which  has  a  cycle  circumscribing 
every  primal  node.  If  we  assign  arc  lengths  to  these  arcs  according 
to  our  convention,  we  will  be  adding  two  more  constraints  to  (9). 
However,  the  previous  optimal  solution  to  (9)  is  also  optimal  for 
this  modification. 

If  we  sum  the  differences  over  all  arcs  of  any  of 

these  dual  cycles  we  will  get  zero  since  each  appears  twice  in 

the  summation,  once  with  a  plus  sign  and  once  with  a  minus  sign. 

Therefore,  X  „  =  V.  -  V.  results  in  the  sum  of  the  flows  entering 
ag  j  i 

and  leaving  a  primal  node  being  equal  to  zero  and  (11)  is  satisfied. 

Finally,  the  optimal  solution  to  (9)  will  result  in 
=  for  all  arcs  on  the  shortest  route  between  nodes  1  and  N 

of  the  dual  and  equation  (10)  will  result  in  the  flows  in  the  arcs 
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of  the  primal  minimum  cut  set  (S,S)  being  M  if  aeS  and  0eS 

ap 

and  L  „  if  aeS  and  8eS.  No  flow  augmenting  paths  will  therefore 
a  p 

exist  between  the  primal  source  and  sink  and  hence  the  existing  flown 
are  optimal.  This  completes  the  proof  cf  Theorem  5. 

The  reader  will  recall  from  Theorem  4  that  we  also  can  get 
the  minimal  flow  if  a  problem  is  feasible  by  changing  the  signs  on  tie 
arc  lengths  of  the  dual  and  finding  the  longest  route  from  the  sink 
back  to  the  source.  The  associated  longest  route  linear  programming 
problem  can  be  stated  as 


maximize 

l  h  Uij ' 

f +1 

(i  -  N), 

N 

1 

subject  to 

£  ^  -  V 

- 

1 0 

(i  -  2,3,... ,N-1)  , 

(12) 

L  -i 

(i-1). 

and  £  0. 

The  corresponding  dual  linear  programming  problem  is 

minimize  W„  -  W. , 

N  1 

subject  to  -  Wj  £  ((i,j)eA).  (13) 

If  we  determine  the  longest  route  from  each  node  to  the 
origin,  we  will  have  an  optimal  feasible  solution  to  (13).  The  value 
of  is  minimal  feasible  flow  through  the  network  (a  negative 

value  means  flow  from  the  sink  back  to  the  source) .  The  optimal 
feasible  solution  to  the  original  primal  minimal  flow  problem  can  then 
be  obtained  from  the  results  of  Theorem  6.  The  proof  of  Theorem  6 
parallels  that  of  Theorem  5. 
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Theorem  6:  ConeideA  a  pnimal  oac  (a, 3)  diAected 
node,  a  to  node.  3  which,  htu  a  llow  x^.  Let  the 
intent  exiting  dual  aA.ce,  be  (i,j)  end  (J,i)  with  the 
length  oh  (i,j)  being  -M  and  the  length  oh  (j,i) 
being  l^.  Jh  xag  taXUhi^ 

Xa3  "  Wj  -  Wi 

ho  A.  eveny  oac  (a,  3)  then  an  optimal  hexuible  solution 
lu u  been  hound  to  the  pAimal  minimal,  blow  pnoblem. 
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5 .  AN  EXAMPLE 


The  results  of  Theorems  1  through  6  will  be  illustrated  by 

the  following  example.  Figure  5a  is  a  network  having  node  1  as  its 

source  and  node  2  as  its  sink.  The  numbers  on  the  arcs  represent  the 

lower  and  upper  bounds  on  arc  flow.  The  dual  construction  is  shown 

by  the  nodes  A,  B  and  C  and  the  dashed  arcs  connecting  them. 

The  resulting  dual  network  for  the  maximal  flow  problem  is  shown  in 

Figure  5b.  The  maximal  flow  through  the  primal  is  4  corresponding 

to  V  -  V  .  The  dual  shortest  route  is  the  chain  of  arcs  in  Figure 
L  A 

5b  having  the  hash  marks.  The  primal  minimum  cut  set  is  therefore 

{(1,2), (2, 3)}. 


4,6  / 


v*»o 


Vb-2 


Vc«4 


Figure  5. 


The  flow  through  arc  (1,2)  of  the  primal  is  6  as  given  by 

the  difference  -  V  since  the  dual  arcs  (B,C)  intersect  it. 

C  B 

The  flow  through  arcs  (2,3)  and  (3,1)  are  both  2  corresponding 

to  the  difference  V.  -  Vn. 

A  B 
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If  is  changed  to  a  value  of  4  then  the  top  (A,B) 

arc  of  Figure  5b  has  a  length  of  -4.  A  cycle  of  negative  length 
will  exist  consisting  of  the  top  (A,B)  arc  and  the  bottom  (B,A) 
arc.  The  problem  is  infeasible  according  to  Theorem  2.  Inspection 
of  Figure  5a  under  this  change  shows  that  an  infeasibility  condition 
has  been  created  by  arc  (3,1)  having  a  lower  bound  exceeding  the 
capacity  of  arc  (2,3). 

Figure  6  is  the  dual  network  for  the  minimal  flow  problem. 


I 


The  minimum  feasible  flow  from  the  source  to  the  sink  of  the  primal 

is  W_  -  W.  *  1.  The  arc  flows  creating  this  are  X10  *  W  -  W  =  4 
L  A  1/  L  d 


and  X 


=  3. 
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6.  MAXIMAL  FLOW  IN  NON -PLANAR  NETWORKS 

We  would  like  to  make  use  of  the  ideas  developed  in  the  preceding 
section  for  all  networks  and  not  just  those  which  are  source-sink  planar. 
The  purpose  of  this  section  is  to  present  a  way  of  approaching  such 

problems . 

Suppose  we  construct  a  two-dimensional  representation  o£  a  non- 
planar  network  such  that  all  arcs  are  straight  lines.  If  we  do  this, 
then  no  two  arcs  will  Intersect  more  than  once  at  points  which  are  not 
nodes.  If  we  designated  these  extra  arc  Intersection  points  as  pseudo 
nodes  we  will  have  created  a  network  which  Is  source-sink  planar  and 
„e  can  construct  its  dual  network  using  the  procedure  and  convention 

of  section  2. 

Unfortunately,  we  have  no  guarantee  that  the  shortest  route 
through  this  dual  network  will  give  the  maximal  feasible  flow  and  that 
the  associated  optimal  feasible  an  'lows  can  be  determined  from 
Theorem  5.  Consider  Figure  7 .  Suppose  arc  (a.B)  and  (y,6)  intersect 

at  pseudo  node  <J> . 


Figure  7. 


23 


If  we  hope  to  use  the  results  of  Theorem  5  to  obtain  the 

optimal  feasible  flows  in  these  arcs,  we  must  require  that  the 

values  V±.  V  ^ ,  Vfc,  and  in  (9)  for  the  associated  dual  nodes 

satisfy 


(14) 


Condition  (14)  is  needed  because  an  optimal  feasible  solution  to  the 
original  primal  maximal  flow  problem  must  have 


Xa*  '  V 
Xy<J>  "  V ’ 

since  arcs  (a.fl)  and  (y,6)  do  not  intersect  in  the  original  primal 
network.  Therefore,  for  Theorem  5  to  be  useful,  condition  (14)  must 
be  added  to  the  constraints  of  (9)  . 

We  note  in  passing  that  (14)  can  be  consolidated  into  only  one 
equation  which  is 

Vj  +  \  "  Vi  "  Vm  =  °*  (15) 

Let  the  node  group  surrounding  a  primal  pseudo  node  <f>  be 
represented  by  (i,j,k,m)^.  Let  0  be  the  set  of  all  pseudo  nodes 
in  a  given  primal  network.  The  incorporation  of  (15)  into  (9)  then 
results  in  the  following  problem: 
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minimize  V„  -  V, 

N  1 

subject  to  V  -  V±  r  l  (U,j)eA),  (16) 

vj  +  \  “  vi  “  vm  "  0  (U.J.MO^i  <M) 

Obviously,  the  optimal  solution  to  (16)  cannot  always  be  deter¬ 
mined  using  an  algorithm  for  finding  the  shortest  route  from  node  1 
to  each  node  i.  The  following  algorithm  is,  therefore,  suggested  as 
a  way  of  solving  (16) . 

Algorithm: 

1.  Compute  the  shortest  route  lengths  from  the  dual  origin  to  all 
nodes.  If,  in  the  process,  a  cycle  of  negative  length  is  detected, 
terminate  the  algorithm.  The  primal  flow  problem  is  infeasible. 

If  no  cycle  of  negative  length  exists  then  compute  the  value  of 
Vj  +  Vk  ”  Vi  "  for  (i,j,k,m)(j)  associated  with  each  pseudo 
node  in  and  go  to  step  2. 

2.  Examine  the  set  of  V.  +  V,  -  V,  -  V  values  obtained  in  step  1 

j  k  i  m 

(a)  If  V,  +  V,  -  V.  -  V  =0  for  all  pseudo  nodes  in  4>  then 

j  k  i  m 

terminate.  The  optimal  solution  to  the  primal  network 
problem  can  be  obtained  by  applying  Theorem  5  to  the  planar 
representation  of  that  network. 

(b)  If  V,  +  V.  -  V.  -  V  #0  for  any  set  (i,j,k,m).  then 

j  k  i  m  <p 

compute  a  set  of  numbers  (i  =  1,2,...,N)  where 

yn  =  vn  and 

Y .  =  max  (Y -l , .) 

1  }  J  iJ 

for  i  =  1,2,...,N-1.  Go  to  step  3. 
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3.  For  each  set  (i,j,k,m)  for  which  V.  +  V.  -  V.  -  V  f  0 

9  j  k  i  m 

compute  the  following  numbers : 


M+>  'V'i’V V' 

t2(*)  - 

If  5,^(4))  <  £2(<f>)  for  one  or  m°re  sets  then  go  to  step  4; 
otherwise  go  to  step  5. 


4.  Determine  that  set  (i.jjk.m)"  corresponding  to 


z2(<p)  -  ^Oj))  *  max{£2((j,)  -  £1(<)))  |  £1(4>)  <  fi_2 (4>) } 


and  change  £^  ,  £  ,  £^  and 

ciated  with  that  set  to 


£  .  of  the  arcs  asso- 
mk 


£ij  =  £km  "  £1((|,)  * 
£ji  =  £mk  =  ’ 


Return  to  step  1. 


5. 


Select  any  set  (i,j,k,m)  for  which  £  (<k)  and 

9  1 

values  have  been  computed  and  change  £^,  £^, 

£  .  of  the  arcs  associated  with  that  set  to 
ink 


i2($) 

£,  and 
km 


lii '  *ta,  ’  h‘*> 
lji  ‘  Eik  ' 


Return  to  step  1. 
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7 .  AN  EXAMPLE 

To  illustrate  the  algorithm  we  will  solve  for  max  Q  in  the 
example  shown  in  Figure  8.  This  network,  while  being  planar,  is  not 
source-sink  planar.  Figure  8  is,  in  fact,  the  two  dimensional  repre¬ 
sentation.  We  see  that  all  arcs  are  straight  lines  and  one  pseudo 
node  has  been  created.  The  numbers  on  the  arcs  are  the  values; 

Ly  =  0  for  all  arcs. 

The  dual  network  is  shown  in  Figure  9.  The  and  values 

are  shown  next  to  their  respective  nodes.  We  see  immediately  that 

\r  +  V  -  V,  -  V  ^  0  so  the  Y.  values  are  needed.  Next  we  obtain 
f  e  d  g  i 

S,^)  =  5  and  (<p)  =  8.  Therefore,  we  change  S,df  and  Si to 

5,  Si,,  and  i  to  -5  and  recompute  the  shortest  route, 
fd  ge 

The  new  set  of  values  are  shown  in  Figure  10.  We  now  have 

V,.  +  V  -  V,  -  V  =0  so  we  are  done.  The  maximal  possible  flow 
f  e  d  g 

through  the  primal  network  is  7  as  indicated  by  the  difference 

V  -  V  .  The  optimal  solution  to  the  maximal  flow  problem  associated 
h  a 

with  the  primal  network  is  provided  by  Figure  11.  The  numbers  on  the 
arcs  represent  The  minimum  cut  set  is  {  (4,5) , (4,2) ,(1,2) , 

(3, 2), (5, 3)}. 


1 
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8.  PROOF  OF  THE  ALGORITHM 

When  V.  +  V,  -  V,  -  V  f  0  for  some  set  (i,j,k,m).  we 
j  k  i  m  *  ’  4> 

must  seek  new  values,  Vl,  V ! ,  V',  V',  which  will  satisfy  the 

i  j  k  m  J 

constraints  of  (16)  .  In  an  attempt  to  avoid  reducing  the  value  of 
VN  -  and,  therefore,  the  value  of  max  Q  we  first  compute  the 
values.  The  differences  Y^  -  Y^  (i  ■  1,2,...,N-1)  are  the 
lengths  of  the  shortest  routes  from  each  of  the  nodes  to  node  N. 
Now 


VN  ' 


yn- 


i’ 


(17) 


since  will  be  less  than  or  equal  to  the  length  of  any 

chain  directed  from  node  i  to  node  N  (equality  holds  only  if 

node  i  is  on  the  shortest  route  from  node  1  to  node  N) .  Because 

Y„  =  V  ,  (17)  reduces  to 

N  N 

Yi  *  V 

Therefore,  if  we  select  Vj  satisfying  the  constraints  of  (16)  such 
that 

Yt  ^  S  (18) 


then  we  are  assured  that  will  not  change. 

Combining  (18)  with  the  inequality  constraints  of  (16)  results 
in  the  following  ranges  for  Vj  -  and  -  V^: 


maxt-Jt^.Yj-V^  2;  Vj  -  V'±  £  minU 

max{-£  .  ,Y  -V,  }  £  V'  -  V.'  £  minU 
mk  m  k  m  k 


irvrV- 


.  ,v  -Y,  }. 
km  m  k 


(19) 
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Let  us  assume  the  ranges  in  (19)  have  some  values  in  common. 

We  next  combine  the  inequalities  of  (19)  with  the  equality  constraints 

of  (16) .  The  result  is  (20) .  Notice  that  l  ,  and  £,.  have  been 

mk  km 


omitted.  This  is  due  to  the  fact  that  Z 


..  m  Z  .  and  £.  *  i l  . 

ij  mi  j  i  mk 


for  every  (i,j,k,m)  as  a  consequence  of  the  dual  construction. 

<P 


■“'-V’VfW s  vj  - v;  ■ v:  -  k  s 


We  realize  that  (20)  reduces  to  the  following  form: 

svj  ‘  vi"  K  V*>- 

when  we  use  the  definitions  of  (4> )  and  from  step  2  of  the 

algorithm. 

Step  4  of  the  algorithm  handles  the  problems  where  the  ranges 
of  (19)  have  values  in  common  for  all  $  which  violate  (15) .  It 
changes  the  arc  lengths  of  arcs  (i,j)  and  (k,m)  to  £^(<J>)  and  the 
arc  lengths  of  (j,i)  and  (m,k)  to  — £  ^  (4> ) .  This  forces  the  condi¬ 


tions 


v»  -  \  -  vj  -  vi  ■ 


V  -  V 
vi  k 


V.  “  Vm  > 

J  m 


while  satisfying  (18).  Note  that  one  node  of  each  pair,  (i,j)  and 
(k,m),  will  have  after  the  new  computations.  Finally, 


V'  -  V' 
N  1 


V  -  V,  so  that  the  value  of  maximal  flow  has  not  been 
N  1 


affected  by  the  new  computations, 


(20) 
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If  the  ranges  of  (19)  have  no  values  in  common  for  some  set 

(i,j,kfm)  then  4.  (4>)  <  £„(4>)  and  the  current  V.T  -  V.  value 
<p  1  l  N  1 

represents  an  infeasible  value  for  max  Q.  Because  must  be 

reduced  by  the  amount 

*2(*>  “  *!<♦)  -  max{i2(4>)  -  A^)  |  A^*)  <  A2(*)} 

before  all  sets  (i,j,k,£)^  can  meet  the  constraints  of  (16),  step  3 

adjusts  the  appropriate  arc  lengths  of  the  set  or  sets  corresponding 

to  <f>  to  At  (4>)  and  -A.  (<J> )  to  provide  this  reduction  in  V„ 

l  i  N 

immediately.  No  further  use  of  step  3  is  necessary  after  this  reduc¬ 
tion.  It  is  possible,  however,  that  adjustments  of  the  step  4  variety 
will  still  be  needed  before  an  optimal  solution  is  reached. 

Changing  the  arc  lengths  of  the  arcs  (i,j),  (k,m)  and  (j,i), 

A  A 

(m,k)  to  A^(<f>)  and  -£^(4>)  values  or  to  £^(<j>)  and  -£^((j>)  values, 

respectively,  results  in  the  optimal  solution  to  the  primal  having  a 

A 

partial  flow  of  £^(ij>)  or  £^($)  in  the  intersected  primal  arc  (y , 6) . 

The  arcs  whose  lengths  have  been  changed  either  remain  as  members 

of  the  tree  of  shortest  routes  from  node  1  to  all  nodes  (the  early 

tree)  or  become  members  in  the  next  round  of  calculations.  Only 

when  all  are  members  of  an  early  tree  will  the  condition  = 

V  -  V.  be  satisfied  if  it  was  not  satisfied  initially, 
m  k 

The  forcing  of  the  reverse  arc  of  a  pair  to  a  value  of  -£^(<{>) 
or  -2,^(($>)  does  not  create  a  primal  infeasibility  (that  is,  a  cycle  of 
negative  length  in  the  dual)  for  the  following  reasons.  Suppose  arc 
(i,j)  was  originally  on  the  early  tree.  The  condition  L  ^  £  M  ^ 
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means  that  the  minimum  possible  value  of  8,  ^  will  be  -M^.  This 
will  occur  only  if  in  the  primal  problem.  If  arc  (i,j) 

was  originally  in  the  early  tree  then  the  shortest  route  from  node  i 
to  node  j  is  over  arc  (i,j).  Any  other  route  in  the  dual  network 
will  have  a  length  at  least  as  great  as  £  =  M^.  ^us*  even 

specifying  A  ■  — does  not  create  a  negative  cycle  with  any  of 
the  routes.  The  worst  that  happens  is  that  a  cycle  of  zero  length 
exists  consisting  of  the  arcs  (i,j)  and  (j,i).  If  arc  (i,j)  is 
now  shortened  to  a  value  of  £^(4)  or  £^(4)  then  it  remains  on  the 
tree  and  £  will  have  a  value  -£^(4)  or  — (4> )  which  is  larger 
than  -M  ^  so  110  negative  cycle  is  created. 

Next,  suppose  arc  (i,j)  was  not  originally  in  the  early  tree. 
There  must,  therefore,  be  some  other  route  between  nodes  i  and  j 
which  is  shorter  (it  might  be  arc  (j,i)).  Reducing  the  length  of 
(i,j)  to  £^(4)  or  £^(4)  causes  it  to  become  an  alternate  member 
of  the  early  tree.  Changing  the  length  of  arc  (j,i)  to  a  length  of 
-£^(<J>)  or  -£.^(4)  creates  a  cycle  of  zero  length  with  arc  (i,j). 

If  the  original  shortest  route  between  nodes  i  and  j  was  not  the 
arc  (j,i)  but  rather  some  other  route  directed  from  node  i  to 
node  j  then  its  length  is  at  least  £^(4)  or  £^(<t>)  or  arc  (i,j) 
would  not  be  the  alternative  shortest  route  after  the  new  computa¬ 

tions.  Therefore,  at  worst  a  cycle  of  zero  length  would  be  created  by 
setting  £^  at  a  value  of  -£^(4)  or  -£^(4).  If,  instead,  the 
original  shortest  route  was  some  route  directed  from  node  j  to  node 
i  then  a  change  in  £  ^  has  no  effect  since  it  is  not  a  member  of  a 
cycle  containing  the  original  shortest  route.  Reduction  of  £  to 
the  value  of  £^(4)  or  £^(4)  will  at  worst  create  a  cycle  of  zero 
length  when  combined  with  this  route. 
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9.  COMMENTS  ABOUT  THE  ALGORITHM 

The  intent  of  the  algorithm  of  section  6  is  to  provide  a  means 
of  solving  dual  networks  for  primal  network  optimal  flows  without  the 
worry  of  whether  the  network  is  source-sink  p’anar.  All  that  is  required 
is  that  the  primal  network  be  drawn  such  that  no  arc  intersects  any  other 
arc  more  than  once.  While  the  algorithm  is  designed  only  to  solve 
maximal  flow  problems,  the  necessary  modifications  for  solving  the 
minimal  flow  problem  are  minor  and  obvious. 

This  algorithm  may  not  be  as  efficient  an  approach  as  the  maximal 
flow  algorithm  of  Ford  and  Fulkerson  even  if  no  pseudo  nodes  are  needed 
because  the  dual  graph  must  be  constructed.  The  existence  of  pseudo 
nodes  could  reduce  efficiency  further  since  several  iterations  of  the 
algorithm  may  be  required  before  an  optimal  flow  is  found  or  the  primal 
network  is  found  to  be  infeasible.  Careful  construction  of  the  two- 
dimensional  representation  of  the  primal  can,  however,  reduce  pseudo 
nodes  to  a  minimum. 

The  algorithm  is  quite  useful  and  efficient  for  sensitivity 
analyses  and  parametric  studies  of  network  flow  as  a  function  of  primal 
arc  capacities  and  for  problems  involving  capacity  expansion  and  reduc¬ 
tion  (see,  for  example,  Doulliez  and  Rao  (1971)  and  McMasters  and  Mustin 
(1970)).  In  such  problems,  the  dual  network  needs  to  be  constructed 
only  once.  The  changes  in  primal  arc  capacities  then  appear  as  arc 
length  changes  in  the  dual  and  a  new  shortest  route  can  be  quickly 
determined  when  necessary  for  each  change. 
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